THE ISENTROPIC EULER SYSTEM ADMITS SOME PLANE WAVE 

SUPERPOSITIONS 
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Abstract. A class of difFerentiable solutions is proved for the isentropic Euler equations in two 
and three space dimensions. The solutions are explicitly given in terms of solutions to inviscid Burgers 
equations, and several directions of propagation. The relative orientation of the directions is critical. 
Within the directional constraints, the Burgers solutions are arbitrary. The several velocities add, 
and the pressures combine nonlinearly. These solutions cannot exist beyond the time when shocks 
develop in any of the Burgers solutions. 
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1. Main result. Consider tlie isentropic Euler equations 

ut + u ■ Vm + p^^\7p — 0, pt + div(/9u) =0, p — kp~' 

We assume 1 < 7 < 3, fc is constant, and set a = 

Theorem 1.1. Letvj be unit vectors mM", d=2 or3, for which the dot products 

Vi ■ Vj =-a, j. (*) 

The number N of such vectors is indicated in the table below. 

Further suppose that fj{s,t) are differentiate solutions to Burgers equation 

/t + (1 + =0, s e M, < t < T 

Define 

N ^ i 

^^fj{x-Vj,t)vj, and p = (-=^fj(x ■ Vj,t)j ° 
j=i ]=i 

Then u and p satisfy the isentropic Euler equations on this time interval and while 

Note that the case 7 = 2 corresponds to the shallow water model and you have 
three vectors Vk coplanar at 120 degrees, while 7 = § corresponds to the monatomic 
gas with four Vk having the symmetry of a regular tetrahedron. 

Proof. We will write out the d = 3 case, and the case d = 2 can be obtained by 
deleting the third component of all vectors. Inspired by the treatment in Lax [3], we 
work with the symmetric hyperbolic form 



qt + Aiq^^+A2qx2+A3q^^^0, q = 

of the Euler equations where q is a 4 x 1 vector consisting of the velocities together 
with w — y^'fp/ p, which is proportional to the sound speed. That gives in the 
isentropic case p — {—^w)i. The two forms of the Euler equations are equivalent 
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for differentiable solutions with p > 0. Here Aj 
identity matrix and 



Ujl + awLj, where / is the 4x4 
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We abbreviate dfj/ds evaluated at {x ■ Vj,t) by fjs- Component i of vector Vj is 
written vji. Also abbreviate dfj/dt{x ■ Vj,t) by fjt, and fj{x ■ Vj,t) by fj. Sums are 
from 1 to A'' = 2, 3, or 4, depending on the number of vectors Vk- 

We will need to know the eigenvectors of linear combinations of the Lj. These 
eigenvectors may be read from the calculation 



h 

A: 

TO 

h k m 



h 




h 


k 


= ± 


k 


TO 




m 


±1 




±1 



ft, 

A; 

TO 

h k m Q 



ho 
ko 
Too 




whenever h + k +m = 1 and hho + kko + mmo = 0. 

Now look for solutions of the form q(x,t) = • Vk,t)zk where constant 

vectors Vk € and Zk € are to be found. Then 



Qt 



E 



{ujl + awLj)q^. = ^ (^fkt + ^{ujl + awLj)fksVkj^Zk 



Now suppose we are looking for eigenvectors J2i{''^kjLj)zk = \kZk- As displayed 



above, we may either choose Zk 



with Afe = ±1, or if Afc = then the first 



three components of Zk must be orthogonal to Vk- 
With any such choices of eigenvectors then 



+ E ^i^'^i = E {f^* + ■{u-Vk + awXkfj Zk 

j k 

= i^fkt + Iks ■ (q ■ )) ^fe = E (•^'^* + E ■ 



Vk 

aXk 



)) 



Zk 



We choose to make the dot products z„ 
system into the equations 



Vk 

aXk 



for fc 7^ to, which decouples the 



fkt + (zk ■ ^fkfk. 



0. 



If Afe = ±1 we have Zk ■ 



Vk 

aXk 



Vk-Vk+ oA| = 1 + o. If Afc = then Zk 







where 



v^ is some vector perpendicular to Vk, and Zk ■ 
of t, as well as v:^ ■ Vk = for k ^ m. 



0, so we need fk independent 



ISENTROPIC EULER SYSTEM 3 

Now we analyze the several cases of dot products and eigenvalues. In the cases 
where some = — 1, we replace Vk by —Vk and /fc(s, t) by — s, t). This effectively 
replaces —1 by +1, and we can assume from now on that all Afc > 0. 

The most important case, and the one stated in the Theorem, is when all eigen- 
values are +1. Let the Vk be N unit vectors with all Vk ■ Vm — —a for k ^ m, and 
all Afc = 1. The N is given in the table. The decoupled equations for the fk are the 
inviscid Burgers [1] equation /fet + (1 + a)fkfks = 0. The solutions are of the form 



Q = EfeLl fk{x ■Vk,t) 



This completes the proof. 



Another possibility is that some eigenvalue is 0. Corresponding to each eigen- 



value you may replace the term fk {x ■ Vk, t) 



by gk{x ■ Vk) 







where is any 



differentiable function, provided that vt is perpendicular to Vk and all the other Vn 
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Fig. 1.1. The table shows the number N of vectors v^. available for various d and 7. 



Remark on the time of existence. Such configurations cannot generally live 
beyond the time when shocks develop in any of the fk ■ For example, suppose a shock 
of speed a develops in /i, and assume 7 = 1.4. The jump condition on density is 
[p\a = [pu] ■ vi or 



But this is not possible. Consider a line segment lying in a plane level set of /s and 
within the shock plane. Along this segment, /2 will in general take a continuous 
range of values, while is constant and /i has different one-sided limits depending 
on the side of the shock plane from which the segment is approached. Since 1/a = 5, 
the jump condition is a polynomial identity in the values of /2. This contradicts the 
fundamental theorem of algebra. 

Preliminary calculations done using clawpack [5], [3] suggest that there is a dis- 
tinction in the appearance of pressure contours in two cases of crossing wave fronts 
shortly after breaking occurs, depending on whether the angles between the fronts 
match equation (*). 
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